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ABSTRACT 

Context. The Planck-HFI all-sky survey from 353 to 857 GHz combined with the IRAS data at 100 pm (3 000 GHz, IRIS version of 
the data) show that the dust properties vary from line of sight to line of sight in the diffuse interstellar medium (ISM) at high Galactic 
latitude (10 19 < Ah < 2.5 x 10 2 ° H/cm 2 , for a sky coverage of ~ 12%). 

Aims. These observations contradict the usual thinking of uniform dust properties, even in the most diffuse areas of the sky. Thus, our 
aim is to explain these variations with changes in the ISM properties and with evolution of the grain properties. 

Methods. Our starting point is the latest core-mantle dust model. This model consists of small aromatic-rich carbon grains, larger 
amorphous carbonaceous grains with an aliphatic-rich core and an aromatic-rich mantle, and amorphous silicates (mixture of olivine 
and pyroxene types) with Fe/FeS nano-inclusions covered by aromatic-rich carbon mantles. We explore whether variations in the 
radiation field or in the gas density distribution in the diffuse ISM could explain the observed variations. The dust properties are also 
varied in terms of their mantle thickness, metallic nano-inclusions, carbon abundance locked in the grains, and size distributions. 
Results. We show that variations in the radiation field intensity and gas density distribution cannot explain variations observed with 
Planck-HFI but that radiation fields harder than the standard ISRF may participate in creating part of the observed variations. We 
further show that variations in the mantle thickness on the grains coupled with changes in their size distributions can reproduce most 
of the observations. We concurrently put a limit on the mantle thickness of the silicates, which should not exceed ~ 10 to 15 nm, and 
find that aromatic-rich mantles are definitely needed for the carbonaceous grain population with a thickness of at least 5 to 7.5 nm. We 
also find that changes in the carbon cosmic abundance included in the grains could explain part of the variations in dust observations. 
Finally, we show that varying the composition of metallic nano-inclusions in the silicates cannot account for the variations, at the 
same time showing that the amount of FeS they contain cannot be > 50% by volume. 

Conclusions. With small variations in the dust properties, we are able to explain most of the variations in the dust emission observed 
by Planck-HFI in the diffuse ISM. We also find that the small realistic changes in the dust properties that we consider almost perfectly 
match the anti-correlation and scatter in the observed f3 - T relation. 


1. Introduction 


Dust is seen everywhere in the Milky Way and it is known 
that dust properties evolv e from the diffuse in t erstellar medium 


(ISM) to denser reg 

ions (Bernard et alJll999t iRidderstad et al. 

2006 

; Schnee et al. 

200| Koehler et al. 20l2 Ysard et al. 

2013 

). Except in a 

handful of studies (see for instance Bot 


et al. 2009 and Veneziani et al. 2010 for studies in emission, 
and Fitzpatrick & Massa 2005, 2007, and 2009 for _studiesMn 
extinction), before Planck-HFI ('Planck Collaboration II 120141) . 
it has generally been assumed that dust properties were the 
same everywhere in the diffuse ISM (see for instance Boulanger 
et al. 1996, for Ah 5 5 x 10 20 H/cm 2 ). However, Planck 
observations seem to contradict this assumption. Indeed, varia¬ 
tions in dust spectral energy distribution (SED) were observed 


by Planck-HFI, which sug 
within the diffuse ISM 

gest that dust properties may vary 

P 

anck Collaboration XXIV 2011; 

Planck Collaboration XVII 

2014; Planck Collaboration Xll 


20l4lPlanck Collaboration XXIXll20l4~ 


Using a pixel-to-pixel modified black-body y 1 fib Planck 
Collaboration XI (2014, hereafter PCXI) modelled the all-sky 
emission from 353 to 3 000 GHz: I v = t Vo B v (T) (v/vof, where 
r Vo is the optical depth at vo = 353 GHz, T is the dust colour 
temperature, and (3 is the spectral index, assumed to be con¬ 
stant from 100 pm to 353 GHz. Besides, PCXI computed the 


dust radiance K = J l v dv, luminosity Lh = 4 jfR /Ah, and 
opacity cr Vo = t Vo /Ah- For the most diffuse areas of the sk>Q 
(10 19 < Ah < 2.5 x 10 2 ° H/cm 2 ), the opacity is found to de¬ 
crease with increasing temperature while the luminosity remains 
almost constant (Lh ~ 3.5 x 10~ 31 W/H). This strongly suggests 
variations in dust properties since variations in the interstellar 
radiation field (ISRF) would lead to an increased luminosity. 
Combining Planck-HFI and Sloan Digital Sky Survey (SDSS) 
data, PCXI derived the dust visible extinction for diffuse ISM at 
high Galactic latitude, E(B - V), to far-IR emission ratio to be 
E(B-V)/‘R ~ 5.4x 10 5 . Additionally, PCXI found that the varia¬ 
tions in T and /3 are globally anti-correlated on the sky. These au¬ 
thors used a two-step fit to limit the effect of instrumental noise 
on the estimate of T and [3 and showed that most of the variations 
measured at high Galactic latitude are truly related to variations 
in the dust emission, even if a significant but subdominant frac¬ 
tion of the observed anti-correlation comes from cosmic infrared 
background anisotropies. The uncertainties on T and (3 due to 
these anisotropies are ~ 0.4 K and ~ 0.1, respectively. 

The goal of the present analysis is to explain the variations 
in the dust SED at high Galactic latitude as presented by PCXI, 


1 These areas correspond to the two areas described as “Low Ah” and 
“Lowest 1%” in the last two lines of Tabs. 3 and 4 in PCXI and have a 
sky coverage of ~ 12%. 
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using a physical dust model dJones et al.ll2013h . These variations 
may be explained by two phenomena. Either the dust environ¬ 
ment is varying (density distribution, radiation field intensity 
and hardness) or dust properties are changing/evolving. In the 
following, we explore these two options: firstly, through build¬ 
ing models of the ISM density distribution and varying the ra¬ 
diation field ; and secondly, through varying the properties of 
the iJones et al.l (120135) dust model. Sect. [5] describes the models 
and tools that we use to describe dust emission and extinction 
in the ISM. Sect. [3] shows the influence of varying the environ¬ 
ment properties on the dust observables, whereas Sect. [4] details 
the influence of dust property evolution. The counterpart of all 
these variations in terms of extinction is also compared to obser¬ 
vations. Concluding remarks are given in Sect. [5] 


2. Dust model and tools 


2. 1. Dust model 


Our start ing point is t he Jones et aD (1201 3h dust model as up¬ 
dated by IKohler et~aJ1 d2014ll . It assumes that dust in the dif¬ 
fuse ISM consists of small grains of aromatic-rich carbon (ra¬ 
dius a < 20 nm) and larger core-mantle grains, for which the 
cores are composed either of amorphous aliphatic-rich carbon or 
of a mixture of amorphous silicates with the normative chemical 
compositions of olivine and pyroxene (with half of the total mass 
of silicate grains in each type). For all types of grains, the man¬ 
tles are made of aromatic-rich carbons with a thickness of 20 nm 
for the carbon cores and 5 nm for the silicate cores. Moreover, 
iron and sulphur are incorporated into the silicate cores in the 
form of metallic nano-inclusions of Fe and FeS. They represent 
10% of the total core volume, of which 30% is FeS and 70% is 
pure Fe. The corresponding elemental abundances and dust-to¬ 
gas mass ratios are listed in Tab. 1. 


In Sect. [4] we deviate from this original composition by 
changing the mantle thickness on the carbon and silicate grains 
and the volume of the FeS inclusions in the silicates. We assume 
that 30%, 50% or 100% of the metallic nano-inclusion volume is 
made of FeS. For the carbonaceous grains, we consider mantle 
thicknesses of 0, 5,7.5,10,20, and 30 nm, and for the silicates of 
0, 5, 7.5, 10, 15, and 20 nm. In all cases, we assume a fixed size 
distribution, meaning that when the mantle thickness increases 
(decreases), compared to the standard fjones et akl (1201 3l) model, 
the size of the core decreases (increases) accordingly. The corre¬ 
sponding elemental abundances and dust-to-gas mass ratios are 
listed in Tab. 1. 


2.2. Methodology 

The optical properties o f the core-mantle grains are derived with 
the BHCOAT routine dBohren & Huffman! 119831) and for the 
variations in the metallic inclusions, we use th e Maxw ell-Garnett 
effective medium theory as described in IKohler et al; (12014i) . 
Once the grain optical properties are known, the corresponding 
emission an d extinction are compute d with the DustF.vQ code 
described in IComniegne et akl (1201 lh . DustEM is also coupled 
to the 3D radiative transfer code CRT (lJuvelall2005b I Ysard et aD 
120121) . This allows us to calculate the dust emission and extinc¬ 
tion for various density distributions and a range of hydrogen 
column densities, Ah, similar to those measured by PCXI. 


2 The DustEM code and associated plotting routines can be down¬ 
loaded at http://www.ias.u-psud.fr/DUSTEM/ 


The variations in the dust observables presented by PCXI are 
based on a pixel-by-pixel modified black-body^- 2 fit. These fits 
were performed for the all-sky IRIS 100 gm band (3 000 GHz, 
Miville-Deschenes & Fagache 2005) and Planck-HFI 857, 545, 
and 353 GHz bands (350, 550, and 850 /mi). To compare the 
results of PCXI to the dust model directly, we fit the modelled 
dust SED in the exact same wa\0, i.e. integration in the IRAS 
and Planck bands, followed by a modified black-body fit. 

The results for the diffuse ISM (10 19 < Ah < 

2.5 x 10 2l) H/cm 2 ) presented in PCXI after computing 7', 
[j and r Vo , can be efficiently summarised with the following few 
parameters: 

- the radiance: K = £ I v dv ; 

- the luminosity: Lh = 4 jfR /Ah ; 

- the opacity: cr Vo = r Vo /A H . 

From these parameters, six figures show the variations and 
relations among dust observables (see for instance Fig. |TJ: 

- Panel a: the /3 - T relation (Fig. 16 in PCXI, for 
which we apply the diffuse ISM mask corresponding to 
10 19 < A H < 2.5 x 10 2 ° H/cm 2 ); 

- Panel b: the Lh vs. T relation (their Fig. 20), in which Lh is 
found to be almost constant; 

- Panel c: the cr Vn vs. T anti-correlation (their Fig 


al so obs erved by Planck Collaboration XVIIl (1201 


iPlanck Collaboration XXIXl (120141) : 

- Panel d: the r Vo to Ah ratio (their Fig. 22) 


20 ), 

and 


Panel e: the E(B 
Panel/: the E(B - 


■ V) to t Vo ratio (their Fig. 23) 
V ) to % ratio (their Fig. 23). 


and 


In the following sections, we compare all our results to the 
observational results of PCXI using these six figures (see for 
instance Fig. []}. In panels a (top left), b (top middle), and c 
(top right), PCXI estimated the hydrogen column density used 
to compute Lh and cr Vo from the Leiden/Argentine/Bonn (LAB) 
survey of Galactic HI for the atomic part (IKalberla et al.l 12005 / 
and from the Planck 12 CO J=0—> 1 map (type 3) for the molec¬ 
ular part, which has a high signal-to-noise ratio and is the 
best sui ted to de tect faint CO regions at high Galactic latitudes 
(iPlanck Collaboration XIIlll2014 ). Consequently, it does not ac¬ 
count for the ionised hydrogen and molecular gas not detected 
in CO that can be present along the lines of sight. As explained 
in PCXI, their estimate of Ah for panels b and c is thus a lower 
limit. For these first three panels (a, b, and c), the observational 
results are the density of points maps, on which we overplot¬ 
ted yellow contours: the central contour encloses 50% of the 
observed pixels and the external contour 75%. Additionally, in 
panels b and c, the median values of the luminosity and opacity, 
in equally spaced bins of temperature, are overplotted as white 
circles with “error bars”, which are their lcr-dispersion. 

In panel d (bottom left), the opacity is plotted as a function 
of Ah, also estimated from the LAB survey and the Planck type 
3 CO map, but smoothed to 30’ resolution. The observational re¬ 
sults are represented by black circles, which are the median val¬ 
ues of the opacity in each bin of column density (equally spaced 
in log), and the black “error bars” are the lcr-dispersion of cr Vo 
in each Ah beam. 

The E(B - V) values used in panels e (bottom middle) and/ 
(bottom righ t) come from the Sloan Digital Sky Survey (SDSS) 
on quasars (iSchneider et al.l 12010 1 These SDSS data cover a 
larger column density range (Ah < 6.25 x 10 21 H/cm 2 ) than 


3 N.B. As shown in Fig. 14 of ljones et alJ 12013 1. the spectral index 
P actually varies across this entire wavelength range. However, fixing 
the /?-value is the best one can do with only four observational bands. 
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Table 1 . Dust model abundances and dust-to-gas mass ratios. The first column lists the aromatic-rich carbon mantle thickness, the 
second column the grain density averaged over the size distribution, and the other columns the elemental abundances (in ppm) of C, 
Si, Mg, O, Fe, and S within the grains. For the carbonaceous grains, we distinguish between carbon in aromatic-rich (mantle) and 
aliphatic-rich (core) forms. The carbons included in the silicates are those of the aromatic-rich mantles. 


Small amorphous carbon grains (a < 100 nm) 


Mantle thickness (nm) p(g/cm 3 ) 

M d /M H 

c 

H 

aromatic 

C 

H 

aliphatic 

C 

H 

total 


0 

1.6 

0.17 x 10~- 


0 

143.5 

143.5 

5 

1.3 

0.14 x 10~ 2 

111.7 


1.7 

113.4 

7.5 

1.3 

0.14 x 10~ 2 

114.6 


0.9 

115.5 

10 

1.3 

0.14 x 10“ 2 

115.4 


0.5 

115.9 

20 

1.3 

0.14 x 10- 2 

117.4 


0.2 

117.6 

30 

1.3 

0.14 x 10~ 2 

118.1 


0.1 

118.2 

Big amorphous carbon grains (100 < a < 200 nm) 

Mantle thickness (nm) 

P (g/cm 3 ) 

M d /M H 

c 

H 

aromatic 

c 

H 

aliphatic 

c 

H 

total 

0 

1.60 

0.64 x 10- 3 


0 

53.1 

53.1 

5 

1.35 

0.54 x 10~ 3 


6.6 

45.0 

51.6 

7.5 

1.33 

0.53 x 10~ 3 


9.3 

41.7 

51.0 

10 

1.32 

0.53 x 10- 3 

11.7 

38.7 

50.4 

20 

1.31 

0.52 x 10- 3 

22.1 

25.9 

48.0 

30 

1.30 

0.52 x 10~ 3 

39.6 


4.4 

44.0 


Amorphous silicate grains (100 < a < 200 nm) 


Mantle thickness (nm) p(g/cm 3 ) M d /M H 


0 

2.95 

0.68 x 10- 2 

44.5 

63.1 

152.1 

29.1 

6.2 

0 

5 

1.60 

0.37 x 10~ 2 

38.5 

54.5 

131.5 

25.2 

5.4 

33.9 

7.5 

1.48 

0.34 x 10~ 2 

36.0 

51.0 

122.9 

23.5 

5.0 

48.1 

10 

1.42 

0.33 x 10- 2 

33.7 

47.7 

115.1 

22.0 

4.7 

61.1 

15 

1.36 

0.32 x 10~ 2 

29.3 

41.6 

100.3 

19.2 

4.1 

85.5 

20 

1.33 

0.31 x 10~ 2 

24.8 

35.1 

84.7 

16.2 

3.5 

111.2 


shown in the four previous figures. These measures were per¬ 
formed along 53 399 separated lines of sight. In panels e and 
/, we only show the PCXI results corresponding to the same 
column density range as in previous figures, 10 19 < /Vn < 
2.5 x 10 20 H/cm 2 , in order to make all of them comparable. The 
observational results are represented by black circles, each one 
being the average of the E(B - V ) values in a r Vo (Fig. e) or R 
(Fig .f) bin. As explained by PCXI, the bin sizes are variable in 
order to keep the same number of samples per bin (N — 1 000). 
The “error bars” are the standard deviation of the E(B - V ) for 
each bin, divided by VN. The solid black lines present the best 
linear fits computed by PCXI for the visible extiction to submm 
optical depth and radiance ratios, respectively. 

Finally, in all panels, the thick coloured lines show our mod¬ 
elling results for a physically reasonable set of dust model pa¬ 
rameters, without attempting to fit the data. 

3. Dust environment 

As suggested by PCXI, one way of explaining the dust observa¬ 
tional results measured in the high Galactic latitude diffuse ISM 
could be to consider variations in the dust environment. Here 
we investigate variations in the gas density distribution in the 
medium, in the ISRF intensity, GiQ. and in the radiation field 
hardness. 

For the density distribution, we consider three cases. These 
include: i) cold neutral medium (CNM): assuming a 2 pc-long 
box filled with a gas with a density varying linearly from 25 
to 70 H/cm 3 ; ii) CNM surrounded by warm neutral medium 

4 Scaling factor for the radiation field integrated between 6 and 
1 3.6 eV. The standard radiation field corresponds to Go = 1 

dDraine & Andersonll 19851) . 


(WNM): considering a 2 pc-long box filled with a 1 H/cm 3 
gas with a denser filament at the centre (1 pc wide and «h = 
100 H/cm 3 ) ; and iii) WNM: assuming three different boxes 
filled with a gas with «h = 1 H/cm 3 and lengths of 20, 30, and 
40 pc. The lengths of the boxes were chosen to match the col¬ 
umn density range measured along the line of sight, 10 19 < N h < 
2.5 x 10 2 ° H/cm 2 , presented in PCXI. These three cloud models 
are th en illuminated by the standard ISRF dDraine & Andersonl 
119851) . for which we vary the intensity from Go = 0.8 to 1.4. As 
seen in Sect. 12.11 for all cases, the dust populations are amor¬ 
phous carbon and silicate grains with aromatic-rich carbon man¬ 
tles 20 and 5 nm thick, respectively. 

The results are presented in Fig. [T] As expected for optically 
thin media, it appears that variations in the gas density distribu¬ 
tion of the medium do not produce any significant changes in 
the dust observables (Figs. Hk b, c, d, e, f, where all the models 
more or less overlap). For instance, the variation in dust temper¬ 
ature does not exceed 0.15 K, when comparing the three cases 
of density distribution. Then, for variations in the ISRF inten¬ 
sity from Go - 0.8 to 1.4, we see that the colour temperature 
varies by ~ 3 K (from ~ 19 to 22 K) but that the spectral in¬ 
dex and the opacity are almost constant (J3 ~ 1.45 to 1.5 and 
cr 353 GHz ~ 8.5 x 10- 27 cm 2 /H). The two highest values of Gq 
(1.3 and 1.4) lead to inconsistent E(B - V)/R ratios. This means 
that even if ISRF intensity variations contribute to changing the 
dust temperature, they cannot explain the observed variations. 
We also investigate the influence of the radiation field hard¬ 
ness on the dust observables. Indeed, hot stars located at lower 
Galactic latitudes and in the Galactic plane may participate in 
the dust heating at high latitudes. These stars have harder radi¬ 
ation fields than the standard ISRF for variable intensities: one 
can expect that when reaching high latitudes, the radiation field 
of these stars would be decreased by a significant factor after 
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Fig. 1 . Influence of the density distribution in the local ISM (silicates with a 5 nm thick mantle and carbons with a 20 nm thick 
mantle). Cold Neutral Medium (CNM) with « H = 25 to 70 H/cm 3 is shown in pink, CNM with = 100 H/cm 3 surrounded by 
Warm Neutral Medium (WNM) with «h = 1 H/cm 3 is plotted in yellow, and WNM with «h = 1 H/cm 3 is displayed in turquoise. 
The three models overlap in the figures. The radiation field is also varied with 0.8 < Go < 1.4. In panels cl, e, and/, seven lines of the 
same colour are plotted. Each of these lines corresponds to a different value of Go- For instance, in panel/, the turquoise lines from 
left to right correspond to Go = 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, and 1.4, respectively. The models falling below the observed E(B - V ) 
to 7? ratio are those illuminated by an ISRF scaled with Go = 1.3 and 1.4. See Sect. l2.2l tor description of the black and white lines 
and symbols. 


going through many insterstellar clouds. To test this, we con¬ 
sider black-bodies with effective temperatures of 10 000, 20 000, 
30 000, and 40 000 K. The black-body intensities are normalised 
to correspond to Go = 1,1.5,2,2.5, and 3. For a given Go, the ab¬ 
sorbed power in the visible/UV range is thus constant, whereas 
the colour changes (i.e. the energy of the average absorbed pho¬ 
ton increases with increasing effective temperature). The results 
are presented in Fig. [2] The model illuminated by a black-body 
at 10 000 K cannot reproduce the observed extinction to radiance 
ratio (orange line in Fig.[2f) but the three other cases fit the data 
with T ~ 19 - 23 K and /3 ~ 1.3-1.45. Variations in the radi¬ 
ation field hardness may thus take part in creating the observed 
variations in the dust observables. 


As can be s een from Fig. EL the standard model of 
iJones et ahl (1201 3l) as updated by Koh ler et al.1 (1201 4t) does not 
fit the median results presented in PCXI (black and white cir¬ 
cles in Fig. Q}. This was to be expected for two reasons. First, 
their model was designed to fit pre-Planck data. Second, these 
authors fitted their model to explain observations of a single line 
of sight in the solar neighbourhood, for which the conditions 
are different from those found at higher Galactic latitude, espe¬ 
cially regarding the column density and the ISRF. The intensity 
of the radiation is indeed fundamental to set the thickness of the 
aromatic-rich carbon mantles on the grains. 


4. Dust properties 

Changes in the intrinsic dust properties are likely to produce 
variations in the dust observables. In this section, we investigate 
the extent of their influence by varying four parameters describ¬ 
ing the grain properties: the thickness of the aromatic-rich car¬ 
bon mantles, the composition of the metallic nano-inclusions in 
the silicate cores, the total abundance of carbon in dust, and the 
grain size distribution. 

The variations in the mantle thickness lead to 42 possible 
combinations between grain populations and the composition 
of the metallic inclusions to 14 combinations (see Sect. 12. 1 b . 
All these combinations are considered when we vary the car¬ 
bon abundance and when we change the size distribution. Thus, 
we explore 280 grain population combinations in total and for 
each, we calculate the emission for seven different ISRF intensi¬ 
ties (i.e. 1 960 emission+extinction models). To keep the figures 
and results as easy to understand as possible, we do not show all 
possible combinations in this section but vary only one property 
at each step. 


4.1. Aromatic-rich carbon mantle thickness 

Grains in t he core-mantl e Jones et all (1201 3l) dust model, and as 
updated bv lKohler et al.ld2014l) . are covered by an aromatic-rich 
carbon mantle, which is 20 nm thick for the carbon grains and 
5 nm thick for the silicate grains. The aromatisation of the outer 
layer of the carbonaceous grains is assumed toresult from photo¬ 
processing by stellar UV-EUV photons dJonesll2012allblf . For the 
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Fig. 2. Influence of the radiation field hardness (silicates with a 5 nm thick mantle and carbons with a 20 nm thick mantle). The 
black-body with an effective temperature of 10000 K is plotted in orange, 20000 K in green, 30000 K in blue, and 40000 K in 
pink. The intensity of the radiation field is also varied with 1 < Go < 3. In panels d, e, and/, as in Fig. [I] each of the seven lines with 
the same colour corresponds to a different value of Go. In panels d and e, the lines are well separated because with the set of model 
parameters tested here, the recovered opacity values are not constant with Go as can be seen in panel c. See Sect. l2.2l for description 
of the black and white lines and symbols. 


silicate grains, this mantle originates from the accretion of small 
carbonaceous particles onto the silicate cores . It was shown that 
the se partic les are likely to be aromatic-rich dBrooke et al .119991: 
iJones et akll2013h . either from the start or by the same aromati- 
sation process as for the carbonaceous particles. The thickness 
of this aromatic-rich mantle may thus vary depending on the 
grain history, age, and environment. Hence we consider a mantle 
thickness varying from 0 to 30 nm for the carbonaceous grains 
and from 0 to 20 nm for the silicates. These thicknesses corre¬ 
spond to abundances of 197 to 273 ppm of carbon and 45 to 
25 ppm of silicon in the grains (the variations in the silicon abu- 
dance come from the fixed size distribution as explained in Sect. 

[IB- 

4.1.1. Mantle thickness on amorphous carbon grains 

The results are presented in Fig. [3] where we assume 5 nm thick 
mantles for the silicate grains. For the carbonaceous grains, we 
consider six different thicknesses: 0, 5, 7.5, 10, 20, and 30 nm. 
For 7.5 to 30 nm thick mantles, the models match the observa¬ 
tions within the lcr-dispersion around the median values. The 20 
and 30 nm cases yield similar results for all the observables (blue 
and pink lines in Fig.[TJ). whereas the 7.5 and 10 nm cases exhibit 
significant differences (green and yellow lines in Fig.Q]). For in¬ 
stance in the case of a 10 nm thick mantle, the spectral index is 
higher by ~ 0.1 (Fig. [3ji, ~ 1.55 or /? ~ 1.65 for a 7.5 nm thick 

mantle), while the opacity at 353 GHz, 7-353 ghz/TVh, is lower by 
~ 20 % and almost perfectly matches the median measurements 
of PCXI for the optical depth and radiance to extinction ratios 
(Figs. [B. d, e). As a consequence, the results presented in the 


following sections are obtained assuming an aromatic-rich car¬ 
bon mantle thickness of 10 nm for the carbonaceous grains. 

We also consider carbonaceous grains without mantles (or¬ 
ange lines in Fig. [3]i. These grains, which are made only of the 
amorphous aliphatic-rich carbon cores, are much warmer and 
their emission peaks at wavelengths shorter than 100 gm. When 
performing the modified black-body fits in the IRIS and Planck- 
HFI bands, the temperature and spectral index of the silicate pop¬ 
ulation are are almost exclusively measured. This leads to higher 
values of/? ~ 1.85 and lower values of T ~ 17.5 - 20.5 K and 
< 7-353 ghz ~ 4.5x 10 - 27 cm 2 /H, which still match the observed dis¬ 
persion but are far from the most often measured values, which 
are the areas surrounded by yellow contours in the figures. For 
an intermediate mantle thickness of 5 nm on top of the carbon 
grains, we find similar model results but closer to the main bulk 
of the observational results (turquoise lines in Fig. [3} /? ~ 1.75, 
cr 353 GHz ~ 5 x 10 - 27 cm 2 /H). This means that relatively thick 
aromatic-rich carbon mantles are definitely needed for the car¬ 
bonaceous grain population to reproduce most of the observa¬ 
tions. 


4.1.2. Mantle thickness on amorphous silicate grains 

The results are shown in Fig. [4] where we assume a 10 nm thick 
aromatic-rich carbon mantle for the carbonaceous grains. Six 
thicknesses are considered for the mantles of the silicate grains: 
0, 5, 7.5, 10, 15, and 20 nm. A 20 nm thick mantle would mean 
that 111 ppm of the cosmic carbon is locked in the silicate grains, 
which means that our models never use more than a total of 
280 ppm of carbon, in agreement with the measurements of cos- 
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Fig. 3. Influence of the aromatic-rich carbon mantle thickness on amorphous carbon grains (for silicate grains with a 5 nm thick 
mantle). Models with no mantle (orange), with mantle thicknesses of 5 nm (turquoise), 7.5 nm (green), 10 nm (yellow), 20 nm 
(blue), and 30 nm (pink) are displayed. In the panels where the orange lines are not visible, it is because they overlay with the 
turquoise lines. The radiation field is also varied with 0.8 < Go < 1.4. See Sect. I2.2l for description of the black and white lines and 
symbols. 



Fig. 4. Influence of aromatic-rich carbon mantle thickness on amorphous silicate grains (for carbon grains with a 10 nm thick 
mantle). Models with no mantle (orange) and mantle thicknesses of 5 nm (turquoise), 7.5 nm (green), 10 nm (yellow), 15 nm (blue), 
and 20 nm (pink) are displayed. In Fig. a, the blue and pink lines are out of the plot with 1.2 < /? < 1.3. The radiation field is also 
varied with 0.8 < Go < 1.4. See Sect. l2.2l for description of the black and white lines and symbols. 
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mic elemental abundances (Parvathi_etaL 2012t lAsnlund et al.l 
l2009HSofia & Meverll200ll:lDwek et aUll 9971) . 


For thicknesses up to 10 nm, the models match the obser¬ 
vations. Thicker mantles (15 and 20 nm) must be rejected as 
their optical depth at 353 GHz is too high by ~ 60 to 90% (blue 
and pink lines in Figs. d) and their /1-values are far from 
the main bulk of the observed dispersion with 1.2 < p < 1.3. 
The extinction to submm optical depth ratios are also too low by 
factors ~ 1.4 and 1.6, for 15 and 20 nm thick mantles, respec¬ 
tively. These results hold when considering combinations with 
carbonaceous grains having 20 and 30 nm thick mantles. For the 
carbonaceous grains without mantle or with a thin 5 nm mantle, 
we can no longer reject the silicates with 15 nm thick mantles. 
This can be understood from Figs. [3je, d. Indeed, the effect of dra¬ 
matically decreasing the thickness of the aromatic-rich layer on 
carbonaceous grains is to decrease the opacity at 353 GHz. This 
is enough to couterbalance the effect of the 15 nm thick mantles 
on top of the silicate grains. However, it seems physically unre¬ 
alistic to assume that the mantles are removed from the carbona¬ 
ceous grains, while, at the same time, increasing their thickness 
on the silicate grains. 


4 . 2 . Composition of the metallic Fe/FeS inclusions 

The iJones et ak l 120131) dust model includes metallic iron nano¬ 
inclusions inside the silicate grains. They are included in the 
form of metal making up 10 % of the grain volume^ and using 
7 0% = 25 ppm of the cosmic iron. Subsequently. lKohler et al.l 
(120141) showed that silicate grains can be a reservoir for inter¬ 
stellar sulphur in the form of Fe/FeS inclusions and that as much 
as 5 ppm of sulphur can be included without major changes in 
the dust emission and extinction (i.e. 30% of the volume of the 
metallic inclusions is FeS). However, the amount of sulphur in¬ 
cluded in the grains is still an open question and, for that reason, 
we investigate the influence of FeS inclusions on the dust observ¬ 
ables. We consider three different volume fractions of FeS in Fe 
inclusions: 30%, 50%, or 100%. This corresponds to 25 ppm of 
Fe and 5 ppm of S, 23 ppm of Fe and 9 ppm of S, and 18 ppm of 
Fe and 18 ppm of S, respectively. 

The results are shown in Fig. [3 The two first cases with 30 
and 50% of the inclusion volume filled with FeS do not show 
any difference (pink and yellow lines, which overlap in Fig. [3). 
However, the case, where 100% of the inclusions are FeS, must 
be excluded as its spectral index is off by ~ 0.1 to 0.2 with 
P ~ 1.25, its optical depth at 353 GHz is too high by ~ 60% 
for Nn > 10 2 ° H/cm 2 (turquoise line in Fig. [33), and its ex¬ 
tinction to optical depth ratio is also too low by a factor ~ 1.5 
(Fig. [3). We conclude that a variable composition of the metal¬ 
lic nano-inclusions cannot account for the variations in the dust 
observables in the diffuse ISM. Conversely, the far-IR/submm 
observations do not provide a strong constraint on the FeS con¬ 
tent (in nano-inclusions) in interstellar dust. 

4.3. Carbon abundance 

In the lJones et al.l (12013l) dust model, about 200 ppm of the cos¬ 
mic carbon is in the grains. The observational estimate of the 
carbon abundance inside the grains varies depending on the ob¬ 
served line of sight. For instance, iParvathi et al.l (1201 2l) found 
that as much as 395 + 61 ppm of the cosmic carbon is in dust 
in the direction of HD 207198. Considering this, we explore the 
effect of adding more carbon in our dust model by making two 
different assessments: first, we add 120 ppm of carbon in the big 


carbonaceous grains (M^/Mh = 1.98 x 1CT 3 ) ; and second, we 
add 60 ppm in both big (M^/Mh = 1.26 x 10 3 ) and small_car- 
bonaceous grains (M^/M H = 2.13 x 10~ 3 ). IParvathi et al1(l2012l) 
also found lines of sight, where the carbon abundance in the dust 
is very low. These observations could be explained by local vari¬ 
ations in the cosmic carbon abundance_or by the destruction of 
carbonaceous grains. Indeed. iBocchio et al.l ( 2014 ) showed that 
in supernova shock waves carbon dust is easily destroyed, lead¬ 
ing to the complete destruction of the smallest grains and to a de¬ 
crease in the amount of big carbon grains. They further show that 
silicate grains are more resilient to shocks than carbon grains. 
Considering this, we explore the effect of including less carbon 
in the grains. We thus remove all the small carbonaceous grains 
and assume that only 12 ppm of carbon are inside the big car¬ 
bonaceous grains (M^/Mu = 0.14 x 1(G 3 ). 

The results are presented in Fig. [ 6 ] for carbon and silicate 
grains with 10 and 5 nm thick aromatic-rich carbon mantles, re¬ 
spectively. For an increase in the carbon abundance in the grains, 
both models match only the optical depth deduced for the lines 
of sight with the lowest column densities ('Nn ~ 6 x 10 1 to 
10 2() H/cm 2 ) but otherwise deviate from the data points (Fig. [ 6 }i). 
Increasing the abundance of carbon in the grains may thus ex¬ 
plain part of the variations measured towards the most diffuse 
areas of the ISM but cannot be the main explanation for the en¬ 
tire range of variations. The results for a decrease in the carbon 
abundance in dust are represented by the green lines in Fig. [3 
For that particular case, we considered radiation field intensi¬ 
ties up to Go = 5. Indeed, if we assume that the carb on grain s 
are destroyed by supernova shock waves (IBocchio et al.ll2014l) . 
the local radiation field is likely to be higher than in the stan¬ 
dard diffuse ISM. We find that decreasing the amount of big car¬ 
bon grains lead to a higher jG — 1.8 and a lower optical depth 
(T 353 ghz ~ 5.5 x 10 27 cm 2 /H, while the luminosity is decreased 
by ~ 10 to 15%. This is thus a likely scenario to explain part of 
the observed variations in the dust properties. 


4 . 4 . Grain size distribution 

The grain size distributions can be altered by proc esses suc h as 
shatte ring and sputtering in the diffuse ISM (IBocchio et al.l2014l 

l20HlJones et d.ll996tlMcKedl989tlSeabll987h . To investigate 

the effect of variations in the grain sizes on the dust observables, 
we test this by dividing the centre radius of the log-normal size 
distributions by a factor of 2 for both silicate and carbonaceous 
grains. 

The results are shown in Fig.[7]for carbon and silicate grains 
with 10 and 5 nm thick aromatic-rich carbon mantles, respec¬ 
tively. The smaller grain populations lead to colour t empe ratures 
about 1 K higher than the standard populations of 1 Jones et al.l 
(12013h . while the spectral index and the opacity at 353 GHz re¬ 
main almost the same (Figs. [7ji and c, p ~ 1.55, T ~ 20 - 23 K, 
and o ~353 ghz ~ 7 x 10 27 cm 2 /H). This allows us to reproduce 
some of the observational results for the diffuse ISM and vari¬ 
ations in the grain size distribution are thus a likely scenario to 
explain at least part of the observed variations. 

4 . 5 . Extinction 

For consistency, we finally check that the extinction curves cor¬ 
responding to our various dust models are qualitatively in agree¬ 
ment wi th the extinction measurements towards th e di ffuse ISM 
made bv lRieke & Lebofskvld 1 985l) . lMathisl (1 1 990l) . and lMcClurel 
d2009h . We are only looking for qualitative agreement as these 
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Fig. 5. Influence of the composition of metallic iron inclusions (for carbonaceous grains with a 10 nm thick aromatic-rich carbon 
mantle and silicates with a 5 nm mantle). We show models, where 30% (yellow), 50% (pink), and 100% (turquoise) of the iron 
inclusions are made of FeS. The first two models almost perfectly overlap in the figures (yellow and pink lines). In Fig. a , the 
turquoise line is out of the plot with (3 ~ 1.25. The radiation field is also varied with 0.8 < Go < 1.4. See Sect. l2.2l for description of 
the black and white lines and symbols. 



Fig. 8. Extinction normalised at 10 gm. The black circles, boxes, 
and the thick black fine are ob serva tional data from iMathisI 
(H990l) . lRieke & Lebofskvl (fT985h . and lMcClurel (i200l . respec¬ 
tively. The solid coloured lin es are for models with the stan¬ 
dard size distributions from [Jones et ahl ( 2013 ). whereas the 
dashed coloured lines show models with smaller sizes as in 
Sect. 14.41 a) The models presented here are for silicate grains 
with a 5 nm thick aromatic-rich carbon mantle and carbonaceous 
grains without mantle (orange), with a 5 nm (turquoise), 7.5 nm 
(green), 10 nm (yellow), 20 nm (blue), and 30 nm (pink) thick 
aromatic-rich carbon mantle, b) The models are for carbona¬ 
ceous grains with a 10 nm thick aromatic-rich carbon mantle 
and silicate grains with no mantle (orange), a 5 nm (turquoise), 
7.5 nm (green), 10 nm (yellow), and 15 nm (blue) thick aromatic 
mantle. 


measurements were done for lines of sight less diffuse than those 
presented in PCXI, with A k ~ 0.1, 0.4, and 0.3, respectively (to 
be compared with our models: 10~ 3 < Ak < 1.5 xl0~ 2 for carbon 
grains with 10 nm thick mantles and silicates with 5 nm mantles 
and 10 19 < IVh < 2.5 x 10 20 H/cm 2 ). The results are presented in 
Fig. [8] for extinction_curves normalised at 10 /mi. The standard 
d ust model of Jones et alJ (1201 3l) reproduces the observations 
of lMathisI (1 1 99()l) . Variations in the aromatic-rich carbon mantle 
thickness of the carbonaceous grains have only a small influence 
on the extinction curve. However, varying the mantle thickness 
on the silicate grains significantly changes the height and width 
of the 9.7 /mi silicate feature. While thicknesses up to 5 nm agree 
with the data points from IMathisI (119908 ). thicker mantles (7.5, 
10, and 15 nm) lead to a broader feat ure, which moves towards 
exp laining the obser vational results of lRieke & Lebofskvl dl985h 
an d lMcClurel d2009l) . 

Fig.[8]also shows the extinction curves for the grain size dis¬ 
tributions presented in Sect. 14.41 The main result is that these 
smaller grains lead to increased extinction at short wavelengths 
(visible/UV) due to an increase in absorption. This explains 
why th ese grains are hotter than those of the standard model of 
iJones et all (1201 3l) as shown in Fig. [7] We conclude that our dust 
model can qualitatively explain the extinction measurements as 
well as most of the variations observed in the dust SED in the 
Galactic diffuse ISM. 


5. Conclusions 

Using the Llones et al . ( 201 30 core-mantle dust model, we inves¬ 
tigated whether the variations in the dust colour temperature, 
spectral index, and opacity observed by PCXI in the most diffuse 
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Fig. 6. Influence of the carbon abundance (for carbonaceous grains with a 10 nm thick aromatic man tle and silicates with a 5 nm 
aromatic mantle). The solid yellow lines show the model with the standard carbon abundance of lJones et al.1 (1201 3l) . the solid pink 
lines with an extra 60 ppm in both the small and big carbonaceous grains, and the solid turquoise lines with an extra 120 ppm in the 
big carbonaceous grains only. The solid green lines show the results for the case, where only 12 ppm of carbon are included in the 
big carbonaceous grains (see Sect. l4.3l for details). The radiation field is varied with 0.8 < Go < 1.4. The dashed green lines present 
the results for the same carbon abundance as the solid green lines but for 1.4 < Go < 5. See Sect. I2.2l for description of the black 
and white lines and symbols. 


areas of the Milky Way (10 19 < Nh < 2.5x 10 20 H/cm 2 ) could be 
explained by varying the dust properties. We also tested for the 
consistency with the dust extinction and the luminosity measure¬ 
ments. Our results are summarised in Fig. [9] where we gather all 
the valid dust models that we present (pink dashed areas). 

We show that the median observations (black and white cir¬ 
cles in Fig. as well as most of the dispersion, can be ex¬ 
plained with little variations in the radiation field and rather 
small assumptions about the variation/evolution of the dust prop¬ 
erties. These include variations in the aromatic-rich carbon man¬ 
tle thickness on both silicate and carbon grains and variations 
in the grain size distributions towards smaller grains than what 
is usually inferred from the observations of the standard dif¬ 
fuse ISM. Smaller grains can for instance explain higher colour 
temperatures for constant radiation field. For the mantle thick¬ 
ness, we found an upper limit of ~ 10 - 15 nm for the silicate 
grains and a lower limit of ~ 5 — 7.5 nm for the carbonaceous 
grains. Thicker mantles lead to lower /3-values for both types 
of grains but to higher submm o pacities . We only considered 
aromatic-rich carbon mantles as Ijones et akl ( 2014 ) showed that 
small carbonaceous particles and aliphatic-rich carbon mantles 
would be aromatised in less than 10 6 yrs in the standard diffuse 
ISM (Ay < 0.7 and Go = 1). Thus, for the high Galactic lat¬ 
itude diffuse ISM, we favour a dust model with aromatic-rich 
carbon mantle thicknesses in the range 7.5 - 20 nm on carbona- 
ceou s grains and 5-10 nm on silicate grains, with the stan¬ 
dard jlon^erar! (12013 ) dust size distributions. Local variations 
in the carbon abundance inside the grains may also take part in 
producing some of the observed variations in the far-IR/submm 


SED. Finally, our dust model can also explain the f3—T relation 
measured by PCXI. We show that the scatter in the /3-values pro¬ 
duced by variations in the dust intrinsic properties are enough to 
explain most of the observed scatter within 18 < T < 23 K and 
1.3 < /3 < 1.8 (Fig. |9ji). Our results thus suggest that the /3 - T 
and <t - T anti-correlations observed in the diffuse ISM could 
simply be due to a scatter in the dust properties. 

In summary, we have shown the capacity of the new dust 
model to explain the observed trends with physically-reasonable 
d ust p roperty variations. We therefore consider the Ijones et akl 
( 2013 ) model to be rather robust. However, along a few lines of 
sight higher temperatures are found (T > 22.5 K), which we 
cannot explain with changes in the dust properties or in the radi¬ 
ation field. These few lines of sight might include HII or highly 
shocked regions, or even grains with peculiar properties due to a 
peculiar history. Detailed analyses of well-known lines of sight 
(radiation field colour and intensity, local density, gas properties, 
carbon depletion, dynamics...) are needed to understand exactly 
which mechanisms produce the observed variations in the dust 
properties, and how much these properties depend on grain his¬ 
tories, ages, and environments. 
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Fig. 7. Influence of the grain size distribution (carbonaceous grains with a 10 nm thick aro matic -rich carbon mantle and silicates 
with a 5 nm mantle). The yellow lines show the model with the standard size distribution of lJones et al.l (12013l ) and the pink lines 
the case where the centre radius of the log-normal size distributions of both carbon and silicate grains have been divided by two. 
The radiation field is also varied with 0.8 < Go < 1.4. See Sect. l2.2l for description of the black and white lines and symbols. 
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Fig. 9. Observational data points are the same as in Figs. 1 to 6. The pink dashed areas show all the values of the dust observables 
that can be produced by the valid dust models presented in Sects. [2 and [4] i.e. all dust population combinations except 100% FeS 
nano-inclusions, silicates with 15 and 20 nm thick aromatic-rich carbon mantles, and aliphatic-rich carbon grains without mantles. 
For models with increased carbon abundances (see Sect. 14.3I >. we included only those with N H < 10 2 ° H/cm 2 . Also plotted are 
models with a decreased carbon abundance. The radiation field is varied with 0.8 < Go < 1.4. We further include models for which 
dust is illuminated by harder radiation fields than the standard ISRF (Sect. [3] for black-bodies at 20 000, 30 000, and 40 000 K and 
Go = 1, 1.5, 2, 2.5, and 3). See Sect. l2.2l for description of the black and white lines and symbols. 


11 





















